Detecting entanglement by symplectic uncertainty relations 



o 
o 

< 



> 
m 
oo 

o 
in 
o 

S3 

-i— > ■ 
G ■ 
cd 

3 : 



Alessio Serafini 

Institute for Mathematical Sciences, Imperial College London, London SW7 2PG, United Kingdom, 
QOLS, Blackett Laboratory, Imperial College London, London SW7 2BW, United Kingdom, and 
Department of Physics & Astronomy, University College London, 
Cower Street, London WC1E 6BT, United Kingdon^ 

A hierarchy of multimode uncertainty relations on the second moments of n pairs of canoni- 
cal operators is derived in terms of quantities invariant under linear canonical (i.e. 'symplectic') 
transformations. Conditions for the separability of multimode continuous variable states are derived 
from the uncertainty relations, generalizing the inequalities obtained in [Phys. Rev. Lett. 96, 110402 
(2006)] to states with some transposed symplectic eigenvalues equal to 1. Finally, to illustrate the 
methodology proposed for the detection of continuous variable entanglement, the separability of 
multimode noisy GHZ-like states is analysed in detail with the presented techniques, deriving a 
necessary and sufficient condition for the separability of such states under an 'even' bipartition of 
the modes. 



I. INTRODUCTION 

The effective detection of entanglement in continuous 
variable systems has become a task of crucial importance, 
notably in view of the rising relevance of such systems 
for quantum information and communication protocols, 
especially in quantum optical settings 0. To this aim, 
several methods have been proposed lately, variously em- 
phasising theoretical and practical aspects of the matter. 

In first seminal contributions to this aim, the positivity 
of the partial transposition has been shown to be neces- 
sary and sufficient for the separability of (1 + n)-mode 
Gaussian states (i.e. of states with Gaussian Wigncr 
function shared by two subsystems of 1 and n modes 
respectively) 0, Q- Furthermore, a general - 'compu- 
tationally efficient', but not analytical - criterion has 
been found for the separability of any Gaussian state 
under any bipartition, based on the iteration of a non- 
linear map |j| . More recently, a full set of conditions has 
been developed to detect genuine multipartite entangle- 
ment in multimode states and a complete theory of 
optimised linear and 'curved' entanglement witnesses for 
second moments (capable of detecting even bound en- 
tangled states) has been developed @. Also, a hierarchy 
of inequalities for higher moments (generally sufficient 
for entanglement) has been obtained allowing to sig- 
nificantly improve the entanglement qualification of non 
Gaussian states as well. 

This paper focuses on the separability of multimode 
continuous variable states under any bipartition. The 
approach of Ref. , introducing a class of (generally) 
sufficient conditions for entanglement on second statisti- 
cal moments, will be recovered and discussed in detail. 
Moreover, such an approach is extended to a more general 
class of inequalities, capable of dealing with 'pathologi- 
cal' instances (entangled states which were not detected 
by the original inequalities derived in Ref. ||). Finally, 
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a relevant class of multimode continuous variable states 
("GHZ-like") will be explicitly considered and analysed 
with the presented methodology, yielding a compact nec- 
essary and sufficient condition for the separability of such 
states of 2n-mode systems under (n + n)-mode biparti- 
tions. R is worth remarking that, as is clearly the case for 
all conditions on second moments alone, the experimen- 
tal test of the inequalities we will derive does not require 
a complete tomography of the state under examination - 
generally a daunting task - but only the measurement of 
the second moments, which might turn out to be delicate 
for multimode systems but seems in the reach of present 
techniques 0. 

The inequalities we will present are based on a class 
of uncertainty relations on second moments, which may 
be compactly expressed in terms of symplectic invariants 
(i.e. of quantities invariant under linear canonical trans- 
formations). Besides retaining a major fundamental in- 
terest per se, uncertainty relations have also lately been 
acknowledged as a main ingredient in the detection and 
qualification of entanglement, whose general scope goes 
well bey ond the scenario of continuous variable systems 
[ToL llll flU Il3j . In fact, for any state g of a bipartite 
quantum system, the positivity of the partially trans- 
posed density matrix g (obtained from g by transposing 
the Hilbert space of only one of the two subsystems) is a 
necessary condition for the state to be separable jlj, [l5j . 
In other words, the violation of the positivity of g is a 
proof of the presence of quantum entanglement in the 
state g, which can in such a case be exploited to various 
quantum informational aims. Now, uncertainty relations 
for quantum observables derive only from the commu- 
tation relations and from the positivity of the generic 
density matrix g. Therefore, any relation derived from 
an uncertainty relation by replacing the state g with the 
partially transposed state g provides a way of testing the 
positivity of g and constitutes thus a sufficient condition 
for the state g to be entangled. In the present instance, 
since the original uncertainty relations come in terms of 
second moments alone, the entanglement conditions are 
cast in terms of second moments as well, and are therefore 
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of great experimental relevance. Furthermore, as it will 
be shown, these conditions turn out to be also necessary 
for the presence of entanglement whenever the positivity 
of the partial transpose is also sufficient for separability, 
namely in (l+n)-mode Gaussian states and (m+n)-mode 
bisymmetric Gaussian states 0, 0, . 

This paper is organised as follows. Canonical systems 
of many modes (like discrete bosonic fields in second 
quantization or motional degrees of freedom of material 
particles in first quantization), notation and uncertainty 
relations are introduced in Sec.|nJ In Sec. lIIII quantities 
invariant under symplectic operations on the field modes 
are constructed as functions of the second moments of 
the field operators. In terms of such invariant quantities, 
simple uncertainty relations for the second moments of 
any n-mode system are derived in Sec. IIVI The par- 
tial transposition of such a relation will promptly lead to 
conditions for entanglement in terms of the second mo- 
ments for any bipartition of the modes in Sec.[V] Finally, 
specific instances of multimode states are considered in 
Sec. I VII and conclusions are drawn in Sec. IVIlj 

II. UNCERTAINTY RELATIONS FOR 
CANONICAL SYSTEMS 

Let us consider a quantum mechanical system de- 
scribed by n pairs of canonically conjugated operators 
{xj,pj}, each of them satisfying the canonical commu- 
tation relations (CCR). As well known, by virtue of the 
Stone- von Neumann theorem the CCR admit, for any 
finite n, only a unique, infinite dimensional represen- 
tation, thus allowing for the occurrence of continuous 
spectra for the canonical operators. Therefore the sys- 
tem in case, whose variables could be motional degrees 
of freedom of particles in first quantization or quadra- 
tures of a bosonic field in second quantization, is com- 
monly referred to as a "continuous variable" (CV) sys- 
tem. Grouping the canonical operators together in the 
vector R = (xi^m, . . . ,x n ,p n ) T allows to concisely ex- 
press the CCR [3 (with h = 2 for future convenience) 
as 

[Rj,Rk]=2in jk , (1) 
where the symplectic form f2 is defined by 

n = ®> , *> = ( _°! J ) • (2) 

A major role in the description of the system is played 
by the transformations which, acting in Heisenberg pic- 
ture on the vector of operators R, preserve the funda- 
mental commutation relations set out in Eq. (|TJ. All the 
conceivable dynamical evolutions of the system have to 
comply with such a prescription. In particular, Eq. 
shows that allowed transformations acting linearly on 
R must preserve the symplectic form Q under congru- 
ence: such linear transformations form the real symplec- 
tic group Sp2n,R (reality being necessary to preserve the 



hcrmiticity of the canonical operators): 

S £ 5L(2n,R) : S e Sfon.R S T QS = Q . 

At the Hilbert space level, these transformations corre- 
spond to unitary operations generated by second order 
polynomials in the canonical operators. This fact, rig- 
orously justified by the construction of the infinite di- 
mensional 'metaplectic' representation of the symplectic 
group (allowed by the Stone- von Neumann theorem), can 
be heuristically understood by considering the Heisen- 
berg equation of motion for R under the action of a 
quadratic Hamiltonian and by recalling Eq. JQ). 

Clearly, any state of a n-mode CV system is described 
by a hermitian, positive, trace-class operator g. Let us 
define the matrix of second moments, or "covariance ma- 
trix" (CM), cr (with entries o-j k ) of the state g as 

a jk ee Tr [(R 3 R k + R k R 3 )g}/2 - Tr [i^Tr [R k g] . (3) 

The CM is a symmetric 2n x 2n matrix. The positivity 
of g and the CCR imply the following semidefinite con- 
straint to be satisfied by any bona fide CM cr 0, [2(j 

cr + in > . (4) 

This well known inequality (whose proof is, for the ease 
of the reader and to provide a self-contained exposition, 
reported in Appendix A) is the only constraint a sym- 
metric 2n x 2n matrix has to satisfy to qualify as the CM 
of a physical state. For future convenience, let us write 
down the CM of an n-mode system in terms of 2 x 2 
sub-matrices ~fj k 

(in ■■■ lln \ 

*= ;••.; ■ ( 5 ) 

\ 1 ' nl ' ' ' Inn / 

Symmetry implies jj k = 7^. Let us also notice that, be- 
cause of the skew-symmetry of ft, Inequality (0} ensures 
the strict positivity of cr |21j : 

cr > . (6) 

For a single mode system, the uncertainty principle Q 
can be more explicitly recast as 

Det a > 1 , (7) 

which corresponds to the well known Robcrtson- 
Schrodingcr uncertainty relation [22I |23| . Inequality {7J) 
together with the condition 10 are equivalent to the un- 
certainty principle for single-mode systems. For a 
two-mode system, the uncertainty principle reads 

Det cr + 1 > , (8) 
where A? ee E^,fe=i Dct 7 jfe - 
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Note that the quantities Det er and , entering into 
Inequalities J7J and @, are invariant under symplectic 
transformations j24|. Clearly, one could provide a com- 
plete set of conditions determining the physicality of the 
CM by just verifying the positivity of the matrix er + ifl, 
which can be checked by considering the determinants of 
the principal submatrices (see, e.g., |25j]1. Still, the in- 
variant nature of the uncertainty principle implies that 
the explicit expression of the uncertainty relation for a 
general n-mode CM er must be possible in terms of sym- 
plectic invariants constructed from the entries of er \ In 
the following section we shall single out such invariants, 
by which an explicit, general expression for a multimode 
uncertainty relation will be subsequently derived. 

III. CONSTRUCTION OF THE SYMPLECTIC 
INVARIANTS 

To begin with let us recall that, because of the pos- 
itivity condition i B Hl , o ne can apply a seminal result 
by Williamson |26l W\. |2S| (originally concerning the 
classification of quadratic Hamiltonians) to the positive 
quadratic form er to infer the following result, which con- 
stitutes a basic tool of symplectic analysis: for any CM 
er there exists a (non-unique) symplectic transformation 
S G Sp2n,TR, such that 

S T (tS = v , 

where 

v = diag(iA, , vj) . 

The quantities {vj} are referred to as symplectic eigen- 
values, while the transformation S is said to perform a 
symplectic diagonalisation of the CM er. The matrix u is 
the so called Williamson normal form of the CM er. The 
uncertainty principle Q| can be succintly recast in terms 
of the symplectic eigenvalues {vj} as 

uj > 1 (9) 

for j = 1, . . . , n. Inequality is fully equivalent to the 
uncertainty relation Q. Its proof is immediate, since 
it derives from the straightforward computation of the 
eigenvalues of the matrix v + iQ, where v is the normal 
form associated to the CM a. 

We can proceed now to the construction of the sym- 
plectic invariants for an n-mode state. The symplectic 
eigenvalues are of course symplectic invariants but, for 
an n-modc system, their analytical expression in terms 
of the second moments turns out to be rather cumber- 
some (when possible at all). As a matter of fact, the 
symplectic eigenvalues can be computed by diagonalising 
the matrix f2er, whose eigenvalues turn out to be {T^j} 
for j = 1, . . . , n. The latter statement is easily proved 
by checking it on the normal form u and by considering 
that 

Civ = QS T aS = S^flcrS for some S G Sp 2 „,Tn . (10) 



Acting by similarity, the transformation S preserves the 
spectrum of f2er which thus coincides with the one of 
ills. Now, a natural choice of symplectic invariants, dic- 
tated by Eq. 1)11)0. is given by the principal minors of the 
matrix Clcr, which are invariant under symplectic trans- 
formations acting by congruence on cr. As we will see 
shortly, a compact and elegant uncertainty principle can 
be expressed in terms of such invariants for any number 
of modes. Let M k(ct) be the principal minor of order 
k of the matrix a j29j. Then, we define the symplectic 
invariants of a n-mode state {A™} for j = 1, . . . , n as 

A™ = M 2j (fla) . (11) 

The principal minors of odd order vanish because of the 
alternate sign in the spectrum of Clcr, leaving us with n 
independent symplectic invariants {A™} (as one should 
have expected, since the number of symplectic eigenval- 
ues is n). The quantities {A™} are also known as "quan- 
tum universal invariants" |30j . We are now interested in 
the expression of the invariants {A™} in terms of the sym- 
plectic eigenvalues {vj}, which can be straightforwardly 
retrieved by considering the normal form v and turns out 
to be the following 

a " = e n . ( i2 ) 

sv- fees™ 

where the sum runs over all the possible j-subsets iS™ of 
the first n natural integers (i.e. over all the possible com- 
binations of j integers smaller or equal than n). In the 
following, A"(^i, . . . ,v n ) will stand for the functions re- 
lating the symplectic eigenvalues to the symplectic invari- 
ants, according to Eq. I|12f> . Clearly, one has A™ = Det er 
while, for two-mode states, it is immediate to verify that 
the invariant A^ coincides with the quantity appearing in 
Inequality (JHJ) (this correspondence is clarified and gener- 
alised in Appendix B, where an alternative way to com- 
pute the invariants {A™} is described). 

The following section is devoted to understanding how 
the uncertainty relation constrains the values allowed for 
the symplectic invariants of a n-mode state, giving rise 
to 'symplectic' uncertainty relations. 

IV. SYMPLECTIC UNCERTAINTY 
RELATIONS 

Let us first introduce some further notation. Let us 
consider an n-mode CV system and define the quantity 

n 

s « - E(- 1 )" + ' A " . ( 13 ) 

3=0 

where we assume Aq = 1. Likewise, E„(^i, . . . ,v n ) will 
stand for the function relating the symplectic eigenvalues 
{vj} to £„, according to Eqs. H12|) and (|13f) . 
Symplectic uncertainty relation: Let cr be the co- 
variance matrix of a n-mode continuous variable state. 
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The symplectic invariant E n; determined according to 
Eqs. and j lMfi . fulfills the inequality 

E„ > . (14) 



Proof. The statement (obviously true for n = 1, for which 
it reduces to v\ > 1) will be proven by induction. Let us 
then assume that, for a generic n, E n _i > 0. One has 

d v ->^n{vi, ■ ■ ■ , V n ) = S n _i(j/i, . . . , Vk-li v k+l} • ■ ■ ) v n) > 

(because of the inductive hypothesis) for any k < n 
and for any value of the {vj}. Therefore, because of the 
bound @, the minimum of E„ is attained for i/j = 1 
Vj < n. In such a case one has 

n 

En = (-1)" ^(-1)^ = 0, 
J=l 

where C™ is the bynomial coefficient. Inequality (|14|> is 
thus established. □ 

Alternately it may be shown, exploiting Eqs. (|12|) and 
(U2J), that 

n 

s„=nK ? -i)- (15) 

J'=l 

Inequality l|14|l . reducing to the well known relations 
and © for, respectively, n = 1 and rt = 2, provides 
a general and elegant way of expressing a necessary un- 
certainty relation constraining the symplectic invariants. 
Rather remarkably, such a relation can be analytically 
checked for any (finite) number of bosonic degrees of free- 
dom. However Eq. l|15|l shows that, actually, Inequality 
(|14f> is only necessary and not sufficient for the full uncer- 
tainty relation (0 to be satisfied as it is not able to detect 
unphysical CMs for which an even number of symplectic 
eigenvalues violates Inequality ©■ 

It is worth remarking that this impossibility is perti- 
nent to the considered uncertainty relation and is not due 
to any fundamental lack of information in the symplectic 
invariants. To better clarify this point, let us first con- 
sider the simple case n = 2: in this instance Inequality 
(|14J) cannot spot the unphysical case Vj < 1 for j = 1, 2. 
However, this situation can be ruled out by verifying the 
additional proviso Deter > 1. Summarizing, the uncer- 
tainty principle is equivalent, for two-mode states, to 
the set of conditions 

A^-A^ + A^>0, A|>1, <t>0. (16) 

(The case Vj < cannot be discriminated by symplectic 
invariants, since they are functions of the squared sym- 
plectic eigenvalues: the additional condition er > has 
thus always to be required.) 

More generally, for n-modc systems, the knowledge of 
the symplectic invariants allows for retrieving the sym- 
plectic eigenvalues, so that the complete uncertainty re- 
lation @ can be in principle recast in terms of sym- 
plectic invariants alone. To show this fact, let v denote 



any one of the symplectic eigenvalues, say v-y, and let 
Al' 1 = A%- 1 (v 2 ,...,v n ). Then one has 

A£ = v 2 Alz{ + A^ 1 for k = 1, . . . , n , (17) 

where A™ -1 = is understood. Once the symplectic 
eigenvalues {A™} are known, the quadratic system H17|) 
for the unknowns v and {A^ 1 } can be solved. Such 
a system of n equations allows for n sets of solutions, 
providing the n symplectic eigenvalues as the n solutions 
for v. This procedure provides an efficient way of deriving 
the symplectic eigenvalues of the CM cr and yields, in 
principle, a recipe to recast the full uncertainty relation 
in terms of symplectic invariants. 

Eq. a 1 of) also allows to elucidate another particular 
facet of the uncertainty relations. Suppose that one sym- 
plectic eigenvalue Vj satisfies Uj = 1. As apparent, in such 
a case the symplectic uncertainty relation (|14J) is satu- 
rated. However, the complete uncertainty relation Q is 
not necessarily fully saturated, as the other symplectic 
eigenvalues could be greater than one. This situation is 
described as partial saturation of the uncertainty rela- 
tion. 'Partially saturating' states have been proven to be 
endowed with peculiar properties. For instance, concern- 
ing two-mode Gaussian states, partially saturating states 
turn out to coincide with the states with minimal entan- 
glement for given purities |3ll l33 | , which in turn prove 
to be relevant in characterizing the multipartite entan- 
glement structure of three-mode Gaussian states |3^ . 

Actually, even for states featuring partial saturation 
a simple necessary uncertainty relation may always be 
found in terms of the quantities {A"}. In fact, let us 
suppose that p < n symplectic eigenvalues Uj are equal 
to 1. In this case, one can iteratively exploit Eq. (|17|) . for 
v = 1, to express all the invariants of lower order {A™ _p } 
for j = 0, . . . , n — p in terms of the original invariants 
{A™}, obtaining 

n—p+l n—p+2 n 

fci=j+i fc 2 =fci+i fcp=fc P -i+i 

. (18) 

where {kj} is a set of p summation indexes. Finally, 
one can recast the uncertainty relation of lower order 
E J (-i)"" P+J A ™~ P > simply as 

n— p n—p-\-l n 

E E ■■• E (-D" +fcp A£ p = 

j=0 k!=j + l fcp=fcp_i + l , , 

n ^ 

E[(-l)" + ^!^A»] >0, 

j=p 

where C J p is the binomial coefficient defined as customary 
by = j — p) lp\}. Eq. (|TT?)) cannot be saturated 
if, as per hypothesis, only p symplectic eigenvalues are 
equal to 1. Therefore, in case of saturation of Inequality 
(|14[1 , one may proceed to test Inequality l|19f) for increas- 
ing values of p, until the saturation disappears: the first 
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nonzero value will determine the legitimacy of the CM. 
Note that, in general, Inequality (|19fl hold for partially 
saturating states, when at least p symplectic eigenvalues 
are equal to 1 (in fact, the inequality has been derived 
under such an assumption). 

In the next section, criteria for the separability and 
entanglement of quantum states will be obtained from 
the uncertainty relation previously derived. 

V. SYMPLECTIC SEPARABILITY CRITERIA 

The positivity of the partially transposed state ("PPT 
criterion" ) is a necessary condition for the separability 
of any bipartite quantum state (i.e. for the possibility of 
creating the state by local operations and classical com- 
munication alone). Conversely, the violation of such pos- 
itivity is a sufficient condition for a quantum state to 
be entangled, in which case quantum correlations are at 
disposal and may be exploited for quantum informational 
tasks. Moreover, as far as the CV systems here addressed 
are concerned, the PPT criterion turns out to be suffi- 
cient as well for the separability of (l + ?i)-mode Gaussian 
states (i.e. of (1 + n)-mode states with Gaussian Wigner 
and characteristic functions) and of bisymmetric (m+n)- 
mode Gaussian states (here and in what follows, we refer 
to a bipartite '(m + n)-mode' CV state as to a state sep- 
arated into a subsystem A of m modes, owned by party 
A, and a subsystem B of n modes, owned by party B). 

These facts come in especially handy for CV systems, 
as the action of partial transposition on covariance matri- 
ces is easily described. Let g be a (to + n)-mode bipartite 
CV state with 2 (to + n)-dimensional CM cr. Then the 
CM cr of the partially transposed state g with respect 
to, say, subsystem A, is obtained by switching the signs 
of the ?7i momenta {pj} belonging to subsystem A. In 
formulae: 

cr = TcrT , with T = 0™ ^ J ^ ) © Ian , (20) 

where t 2n stands for the 2n-dimensional identity matrix. 
Now, in analogy with Inequality derived from the 
positivity of the density matrix g, a (generally) sufficient 
condition for separability derived by the PPT criterion is 
given by @ 

a + iQ, > (21) 

or, in terms of the symplectic eigenvalues {£,•} of the 
partially transposed CM cr (whose normal form will be 
henceforth denoted by £>), as 

Vj>\. (22) 

The results of the previous section allow one to recast 
such separability criteria for (m+n)-mode states in terms 
of partially transposed symplectic invariants {A™ + ™}, 
defined by Af +n = My (Oct). 



This simple accessory result will be precious in what 
follows: 

Little lemma. Let cr be the physical CM of a state of 
a (m + n)-mode CV system, with m < n. Let cr be the 
partial transposition of cr with respect to any of the two 
subsystems. Then, at most to symplectic eigenvalues Vj 
of cr can violate Inequality 122\) . 

Proof. Suppose to perform the transposition in the to- 
mode subsystem: cr = T m erT mi where T m is the partial 
transposition matrix in the TO-mode subsystem (equal 
to the matrix T of Eq. Let V(a) be the dimen- 

sion of the subspace upon which the generic matrix a 
is negative definite. Since T m reduces to the identity 
on a (2n + m)-dimensional subspace, Inequality H21(l re- 
duces to the (definitely satisfied) Inequality J3J on such 
a subspace, thus implying T>(cr + iQ) < m. One has 
then D(v + iQ) = V(cr + iQ) < to, where the equal- 
ity holds because the signature is preserved under con- 
gruence transformations ('Sylvester's inertia law') and 
O+iQ = S T (cr+iQ)S for some S € Sp2( m +n),u- Straight- 
forward computation shows that the eigenvalues ofO+iQ 
are given by {pj T 1}, thus proving the result as the {vj} 
have to be positive (cr = T m crT m > because cr > 0). 
The choice of the transposed subsystem is not relevant, 
since (T n T m )(cr + ifl)(T n T m ) = cr' — iQ (where T n and 
&' are, respectively, the partial transposition matrix and 
the partially transposed CM with respect to the n-mode 
subsystem), and cr — iQ > is equivalent to cr +iQ > 0. 
□ 

In analogy with Eq. I|13fl . let us now define 

m+n 

%n+n - £ (_i r +»+.?A™+™ . (23) 
3=0 

The inequality 

Sm+n > , (24) 

being necessary for Inequality l|22|) to be satisfied, is a 
necessary condition for separability under (to + n)-mode 
bipartitions and is thus a sufficient condition to detect 
entanglement in such a multimode system, irrespective 
of the nature of the state under examination. 

Let us remark that, while for Gaussian states Inequal- 
ity (|22|l (which ultimately reduces to a condition on the 
second moments) is necessary and sufficient for the pos- 
itivity of the partial transpose, this is not the case for 
non-Gaussian states. For such states conditions resort- 
ing to higher moments lead to a sharper detection of 
negative partial transposition, and thus of entanglement 


Even for Gaussian states, Inequality (|22|l cannot de- 
tect the negativity of the partial transpose if an even 
number of symplectic eigenvalues violate condition 122|l . 
However, let us focus on Gaussian states under (1 + n)- 
mode bipartitions (for which the PPT criterion is neces- 
sary and sufficient for separability 0). Because of the 
previous lemma, for these states at most one partially 
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transposed symplcctic eigenvalue can violate Inequality 
l)22p. Inequality (|24[l is then capable of detecting such a 
violation. 

The same argument applies to 'bisymmetric' Gaussian 
states, defined as the (TO+n)-mode Gaussian states which 
are invariant under mode permutations internal to the m- 
mode and n-mode subsystems. A bisymmetric Gaussian 
state with CM er can be reduced, by local symplcctic 
operations (on the m-mode and n-mode subsytcms), to 
the tensor product of a two-mode Gaussian state and 
of uncorrelated thermal states with global CM er 2 : 
er 2 = SjcrSi for some Si € Sp2m,R® Sp2n.R- The lemma 
above can be applied to obtain T>(TcrT + ifl) < 1 from 
which, observing that 

(TSjT)(TcrT + in)(TSiT) = Ta 2 T + ifl , 

one infers that at most one partially transposed symplcc- 
tic eigenvalue of the CM er can violate Inequality l(22j) . 
Notice that the locality of the operation Si is crucial in es- 
tablishing this result, since it implies (TSjT)tt(TSiT) = 
n (which would not generally hold for a Si with nondi- 
agonal terms relating the m-mode to the n-mode subsys- 
tem) . 

Inequality (|24|1 is thus necessary and sufficient for the 
separability of all (l + n)-mode and bisymmetric (m + n)- 
mode Gaussian states, except for the set of 'null measure' 
for which any i)j is identical to 1 (that could be entangled 
but have £ m+ „ = 0). However, as in the case of the 
uncertainty relation, a necessary and sufficient condition 
for separability in terms of the quantities {A.,} may be 
found for such states as well. If p symplcctic eigenvalues 
Vj (for p < n) are equal to 1, such a relation reads, in 
analogy with Inequality (|19|L 



m+n 



J2 [{-l) m+n+1 p\C^A 



m+n 



> o 



(25) 



3=P 



Eq. 125(1 cannot be saturated if, as per hypothesis, only p 
eigenvalues were equal to 1. Therefore, in case of satura- 
tion of Inequality (|24f) . one may proceed to test Inequality 
(|23|l for increasing values of p. until the saturation disap- 
pears: the first nonzero value will determine the positiv- 
ity of the partial transposition of the state at the level of 
second moments. Notice that the case for which Vj = 1 
for some j is indeed very relevant in designing and de- 
tecting entangled resources for quantum information as, 
typically, it corresponds to pure CV states. For p = 0, In- 
equality H25f) correctly reduces to Inequality l|24|) . which 
has thus been generalised. Summing up, Inequality \2J$ 
- together with Inequalities (|25|) in case of partial satura- 
tion - are necessary and sufficient for the separability of 
(l + n)-mode and of bisymmetric (m + n)-mode Gaussian 
states. 

The simple condition (|24|) on the second moments, 
which can be promptly analytically verified, may be very 
helpful in the detection of interesting CV entangled states 
(in the same spirit as in Refs. H,E1)- This aspect is espe- 
cially relevant in view of the recent experimental devel- 



opments in the implementation of multipartite CV quan- 
tum information protocols, usually relying on symmetric 
resources (discriminated by the previous condition) |34j . 

The next section will better illustrate the efficacy of 
Inequality l|24|) for the detection of CV entanglement, 
by addressing the separability of a relevant example of 
multimode states. 



VI. EXAMPLE: NOISY GHZ-LIKE STATES 

The CV "GHZ- type" states j3^ are a class of fully sym- 
metric multimode Gaussian states introduced in Ref. [sf3 | 
as the prototypical resource for the implementation of a 
CV teleportation network. Experimentally, they can be 
generated by inserting squeezed vacua into an array of 
beam splitters (for a detailed description of the generat- 
ing scheme, see and |13). Moreover, they turn out 
to be the symmetric Gaussian states maximising both 
the couplewise (between any pair of modes) and the gen- 
uine multipartite entanglement [s^L l38j | and their mul- 
tipartite entaglement proves to be especially resilient to 
decoherence j33|. We will study here the separability 
of GHZ-type states determining the partially transposed 
symplcctic invariants and employing Inequality (|24|l . 

To fix ideas, we will consider a GHZ-type state with 
an even number of total modes In under the n + n bi- 
partition. The CM er of such a state, generated from 
impinging squeezed vacua with squeezing parameter r, is 
described by Eq. (J3J) with 2x2 blocks given by 



1 < j < 2n , 
(26) 




J jk = a 



sinh(2r) 
n 

o 





sinh(2r) 



Realistically, we have assumed a thermal noise with mean 
photon number (q — 1) > to affect the creation of the 
states (this amounts to multiplying the CM by q). 

After some basic algebra (see also Appendix B), 
the partially transposed minors AJ +n can be explicitly 
worked out and found to be 

A n+n = q 2j [(C 2„ _ g n } + g n cosh(4r) ] _ (27) 

where the coefficients g™ are determined by the following 
sums 

9? = fe^^^^^^^^Vfe^ 1 )^*" 1 ^"^ 2 ^" 



(28) 

and satisfy a recursive relation which will shortly be use- 
ful 

g? = 2(C£ 1 -g?_ 1 )- ! #_ 3 for2<i<2n. (29) 
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Also, one has g$ — and = 2. Inequality l|24[l then 
reads 

2n 

Y / (-iy +2n q^ [(C] n - gf) + g? cosh(4r)] > . (30) 

3=0 

Now, the polynomial in q in the LHS of Inequality i|30|) 
admits An — 4 degenerate roots for g = =pl and 4 more 
roots for g = =p e^ 2r [as can be shown by employing 
Eq. H29|) ]. Because the polynomial is obviously diverg- 
ing for infinite q (the leading order in q has a positive 
sign) and accounting for the physical condition q > 1, we 
find that Inequality (|30|l is equivalent to 

q > c 2r . (31) 

Such a simple criterion, derived from Inequality Ij24[l. is 
necessary and sufficient for the separability of 2n-mode 
noisy GHZ-like states under (n + n)-mode bipartitions, 
regardless of the total number of modes (let us recall that, 
in the absence of noise, such states are already known to 
be always inseparable under any bipartition). Note also 
that the physical significance of the derived condition is 
immediately evident: in order to maintain the entangle- 
ment in the final state, the thermal noise has to be coun- 
terbalanced with a corresponding level of squeezing in 
the initially uncorrelated input modes. 

Notably, an analogous analysis - based on Inequality 
(|24|l - may be applied for any number of modes under 
any bipartition and for more general states, getting more 
and more useful as the number of modes increases. 



VII. DISCUSSION AND OUTLOOK 

The separability criteria here presented, compactly 
cast in terms of symplectic invariants, have been proven 
to be necessary and sufficient for two relevant classes of 
Gaussian states and have led to remarkably simple con- 
ditions for the separability of GHZ-like states. Besides 
this immediate usefulness, the presented analysis, deriv- 
ing from a recasting of the uncertainty relations in terms 
of symplectic invariants, touches the very core of the 
symplectic structure, encompassing the full description 
of Gaussian states of light and discrete bosonic systems 
in general. More specifically, it further clarified the spe- 
cific constraints imposed by quantum mechanics on the 
second moments of canonical operators, in a naturally 
canonically invariant framework. 

This fundamental analysis might yield further inter- 
esting results concerning the entanglement characteriza- 
tion of CV states. In particular, the parametrisation of 
Gaussian states through symplectic invariants has pro- 
vided remarkable insight into the entanglement prop- 
erties of two-mode states 0, |33| and could be, em- 
ploying the techniques here presented, carried over to 
the analysis of multipartite continuous variable entan- 
glement, which has been lately drawing considerable at- 
tention mSmiMmElElil mostly in view 



of the remarkable experimental perspectives uprising in 
quantum optical systems. 
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APPENDIX A: PROOF OF INEQUALITY © 

Let us consider the operator y = (R — Tr [gR]) T Y, 
where Y G <C 2n is an arbitrary complex vector. The 
positivity of g implies Tr [gy^y] > 0, from which 

Tr[ey f y] =Y^tY = Y ] (a + iQ)Y >0, (Al) 

where the matrix t has entries Tjk = Tr [gRjRk] — 
Tr [gRj] Tr [gRk] and the CCR {Q have been employed. 
Due to the arbitrarity of Y G <C 2n , Inequality IjAlf) cor- 
responds to Inequality ipfl. which is thus proven. 

Notice how the uncertainty relation has been derived 
assuming solely the CCR and the positivity of the 
density matrix g. 

APPENDIX B: COMPUTATION OF THE 
SYMPLECTIC INVARIANTS 

An alternative way of computing the symplectic in- 
variants defined by Eq. (|12|) is here outlined, providing 
some deeper insight into the symplectic architecture of 
covariances of bosonic systems. 

Let m = {mi, . . . , m k } stand for a ordered A:-subset 
of natural integers smaller or equal than n, such that 
rrij G N and n > rtij > rrij-i for j = 1, . . . , k, and let Af% 
be the set of all such fc-subscts. Now, for l,fh G A/J, let 
us define the 2k x 2k submatrix <jj _ of the CM er as 

^ llimi ' ' ' TiiTOjt \ 

; ■■. ; • (bi) 

where the 2x2 submatrices ~fj k are defined as in Eq. JSJl. 
The symplectic invariants AJJ are then given by 

a»= J2 Dct -L' ( B2 ) 

where the sum runs over all the possible ordered k- 
subset. This definition amounts to considering the sums 
of the determinants of all the possible (also non princi- 
pal) 2k x 2k submatrices of er obtained by selecting all 
the combinations of 2 x 2 blocks (each block describing 
one mode or the correlations between a pair of modes). 
For instance, for any n, one has A™ = J^™ k=1 Dct ~fjk, 
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that is just the sum of the determinants of the 2x2 
blocks themselves (this is the multimode generalisation 
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